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Abstract 

Maxwell and Dirac fields in Friedmann-Robertson- Walker spacetime is in- 
vestigated using the Newman-Penrose method. The variables are all separable, 
with the angular dependence given by the spin-weighted spherical harmonics. 
All the radial parts reduce to the barrier penetration problem, with mostly 
repulsive potentials representing the centrifugal energies. Both the helicity 
states of the photon field see the same potential, but that of the Dirac field 
see different ones; one component even sees attractive potential in the open 
universe. The massless fields have the usual exponential time dependencies; 
that of the massive Dirac field is coupled to the evolution of the cosmic scale 
factor a. The case of the radiation filled flat universe is solved in terms of the 
Whittaker function. A formal series solution, valid in any FRW universe, is 
also presented. The energy density of the Maxwell field is explicitly shown to 
scale as a~^. The co-moving particle number density of the massless Dirac 
field is found to be conserved, but that of the massive one is not. Particles 
flow out of certain regions, and into others, creating regions that are depleted 
of certain linear and angular momenta states, and others with excess. Such 
current of charged particles would constitute an electric current that could 
generate a cosmic magnetic field. In contrast, the energy density of these 
massive particles still scales as a~^. 
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1 Introduction 



The Newman-Penrose|lj formalism of projecting vectors, tensors and spinors onto 
a set of null tetrad bases has proved to be an immensely useful tool with which to 
investigate the properties of quantum fields in curved spacetime. It has been used 
successfully in various black hole geometries [3 El SI- The method has also been 
used more recently[S] to further study the behaviour of Dirac particles in different 
geometries. As the geometry of our Universe is of the Friedmann- Robertson- Walker 
type, it is important to study how electrodynamics is altered by the real expanding 
Universe in contrast to that in Minkowskian spacetime. The behaviour of the quan- 
tum fields in FRW spacetime must be indicative of the exact nature of the geometry 
of our Universe, particularly whether it is flat, closed or open. Properties of matter 
fields in general, and the massive Dirac fields in particular, must have profound con- 
sequences on the structure formation process. Propagation of electromagnetic waves 
in FRW spacetime was studied by Haghighipour[2i (see references therein), among 
other authors. The Dirac field was studied in the NP formalism by Zeccapj and the 
work was continued further by Sharif [H]. This paper is another attempt to use the 
NP method to study the behaviour of Maxwell and Dirac fields in the Universe. A 
slightly different null-tetrad bases are chosen in this paper, resulting in equations 
that take on the familiar barrier penetration forms without much ado. The Maxwell 
equations for free electromagnetic waves are exactly solvable for all cases in terms 
of well known special functions. The spatial part of the Dirac equations are also 
exactly solvable; the time evolution however is coupled to the evolution of the scale 
factor, can be solved in terms of known functions only for the case of a radiation 
dominated universe. 

The FRW line element is usually written in the comoving coordinates as0 

ds" = de - a\t) (y^^ + r'^ {d0' + sin^ ^#')) , (1) 

where a(t) is the scale factor; the curvature constant K can be scaled to the values 
0, +1 or —1 to describe the flat, closed or open universe respectively. In terms of 
the conformal time V = J ^ change of variable r' = r/ (1 + Kr"^ /A) = r/f{r), 

the line element Eq. ^ becomesfTfl] 



ds^ = a^(?7) 



d^^ - -±- (dr' + ^ gi^2 2^^ 

P{r) 



(2) 



Following Chandrasekhar's guidelines for the black hole|TT|, we can introduce the 
null tetrad = (l,-j,0,oj ,n^ = Y (l> J'O'O) '^^^ = ^ (0, 0, -1, -i sin^), and 
the complex conjugate m* . The non vanishing spin-coefficients are /? = —a = 

l^r^^'^^P = -i{ff-^-lf)^^' = iifj+^-lf)^ ^^d7 = -^,where the 
prime and the overdot denote derivatives with respect to r and t] respectively. 

Besides the vanishing of the spin-coefficient k because the / vector forms a congru- 
ence of null geodesies, and of vr telling us that the tetrad bases remain unchanged as 
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they are parallely propagated along / as in Ref.^, the present choice of the bases are 
affinely parameterized so that e also vanishes. The directional derivatives can then 
be written down as -D = = ^D_,A = n^^d^ = —^D^, 6 = m'^dn = -i—L_, 

and 5* = 171*^^3^ = ^L+,where D± = ff^T |^and L± = ^ ^ In terms of 

what may be called the "tortoise" coordinate defined by 



, dr 2 VKr I ^\ , , ^ ~ ° 

= I — = ^=arctan — - — = ^ 2 arctan(r/2), K = 1 , (3) 

2arctanh(r/2), K = -1 



D± = '£r T The tortoise coordinate shown in Fig. [T] reveals that, near the 
origin all the r*'s coincide with each other; also, from Eq. Q it is seen that as 
r ^ 0,r=K ^ r for both K = ±1; thus, at small distances both the closed and open 
universes look basically fiat. 

The Maxwellian wave equations are solved in the next Section, and the solutions 
are used to investigate the behaviour of the energy-momentum in the expanding 
Universe. Section 3 is devoted to Dirac equations; the spatial parts of the four 
components are easily solved, and so is the time dependence of the massless field; 
the time dependence of the massive field is more complicated as it is coupled to the 
evolution of a; nonetheless, the case of the radiation filled fiat universe is found to 
have a solution in terms of the Whittaker function. A formal series solution that 
can be used in any FRW universe is also written down. The behaviour of the Dirac 
current and energy-momentum is then investigated. The final Section makes some 
concluding remarks. 



2 Maxwell Field 

In the NP formalism, the six independent components of the antisymmetric Maxwell 
tensor F^i, are replaced by the three complex scalarspjj $o = F^J^fn'^, $i = 
^Ffj^y {l^^rf + rrf^rrf) , $2 = F^yirrf^^n'' . Then the Maxwell equations become 

(D-2p)$i = (5* + TT - 2a) $0 - K$2, 

(5-2r)$i = (A + /i - 27) $0 - a$2, 

{D-p + 2e)^2 = (5* + 27r) $1 - A$o and 

((5-r + 2/3)$2 = (A + 2/i) $1 - z/$o. (4) 

With the spin-coefficients and the directional derivatives from Sec. 1, and using 
the fact that {l/f)D± [a"(//r)'"$] = a"(//r)™ (^^D± T ^ + ^^^) , Eqs. 
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become 



D_ 



(L_ + cot I 



ar 
1 
ar 
1 
2ar 

T 

2ar 

T 



(L+ + cot e) 



ra 



ra 

T 



and 



(5) 



The first of the two pairs of Eqs. (jS)) can be used together to eliminate $i and be 
combined into a second order equation for $o only, and similarly from the second 
pair for the equation for $2 only. Writing = ^ora^/f and $_ = 2$2or//,the 
resulting two equations take the form 



P 



+ [L± + cot^] $± = 0. 



It is clear that the Maxwell scalars in Eq. 
Then the angular parts become 



are separable into $_ 



(6) 

(r,77)r±i(^,(^). 



' 1 ^ • ^ ^ _|_ 2i cos 6 d 
^sm 9 89 89 sin^ 9 8ip 



+ 



8' 



sin^ 9 8lp^ 



cot^ 9-1\Y- 



^±1 



-\Y- 



±1, 



(7) 



where A is the separation constant. Same eigenvalue A is used for both the solutions 
Kj-ibecause the equation for y+i goes over to that for Y_i with the change of variable 
9 ^ Tx — 9; the two just represent the two helicity states of the photon. The solutions 
are the spin weighted spherical harmonics of Goldberg, et al. jI2| with A = /(/ + 1); 
/ is the total angular momentum of the field including its spin, so in this case for 
photons, / = 1, 2, 3, .... It is obvious that the azimuthal part is just e*™^'^. The radial 
and time parts then become 



:D 



(8) 



The relative normalization between the two helicity states is set by the easily veri- 
fiable Teukolsky type relations 



A0, 



(9) 



Substituting Z = j^D±(l)± for the dependent variable in Eq. ((HI), operating on it 
with D±, and writing = D±D^ = ^ — we end up with the equation of 
motionZ^^Z = VZ, with the potential V = X{f /rY independent of 17. Consequently, 
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the time dependence e"*'^'' can be factored out. For the stipulated time and az- 
imuthal dependence, D± = ^ioj and L± = ^ ± Then the radial parts of 
0±, which may be written as R± in Eq. ((HI), satisfy complex conjugate equations. 
Finally, we are left with the equation for the radial dependence of Z as. 



R = VR, (10) 



which has taken the familiar form of the one dimensional quantum mechanical barrier 
penetration problem with the potentials for the three cases given by 

r /(/ + l)/r2, K = 
V=l l{l + l)/sm^n, K=l . (11) 
[ /(/ + l)/sinhV,, K = -1 

Eq. (fTn| is basically the same as that derived in Ref.jHI. The potentials barriers 
represented by Eq. (fTTll are just the rotational kinetic energies. They become thicker 
for larger /. The barriers go down as 1/r^ in the flat universe, fall off even more 
rapidly in the open universe, and are potential wells between < r^, < vr that rise 
to infinity towards both the ends in the closed universe (Fig. If r^, is replaced by 
ir^ and u by iuo in the potential for i^' = 1 in Eq. (fTTll we end up with the one for 
K = —1. Both the spin states of the photon satisfy the same equation with the same 
real potential. The inversion relation is provided by the normalization condition Eq. 
© as 

Ai?± = D^R. (12) 



2.1 Solutions 

All the three cases of Eq. (fTn|) can fortunately be solved. For K = 1, f{r) = 1 and 
r^, = r. Then Eq. (fTTHl has the general solution jOl 

R = r [Aoji (tur) + BoVi {cur)] , (13) 

where ji and yi are the respective spherical Bessel functions. As / > 1, only the 
j's are regular at the origin. In the asymptotic region, as 2ji = h[^^ + hi — > 
[{—iY'^^e^'^^ + (i)'+^e~*'^''] /u, we can use the inversion relation Eq. (fT2|l to see that 
the incoming waves are dominated by i?+and the outgoing by 

In the closed universe, it is most convenient to write the solutions in terms of 
the Gegenbauer functions C" as 

R = Ai (sinr,)'+'C^t,_i (cosr,) + (sinr,)"' Cjj, (cosr,) . (14) 

Again, the first one is regular within the allowed range of r^,, and u > I + 1. These 
em waves form standing waves. 

We may write the solution in the open universe as 

R = A_ie''^'-*2i^i (^-^,^ + 1,1-^^,^(1 -cothr,)^ 

+S_ie''""* 2F1 (-1, l+l,l + iuj, i(l - cothn) j , (15) 
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consisting of the outgoing and incoming parts at r* ^ cxd. Alternatively, we may 
make the replacements n ir^, and a; — > iu; in Eq. (fTljl to study the behaviour 
near the origin. Once $o and $2 are known, Eq. ^ can be used to determine $1 
and complete the solutions for the electromagnetic field. 



2.2 Energy Momentum 

The energy-momentum tensor of the Maxwell field in terms of the $'s is given by 

4:TrTap = |$o|^^a% + |'^'2|^U/3 + 2 |$i|^ +m(„m^)} 

-4$*$in(„m^) - 4<l>t$2/(„m^) + 2<l>*<l>2m(«m^) + cc, (16) 

where the brackets around the indices of the null-tetrad denote symmetrization and 
cc is the complex conjugates of preceeding terms; T^/j is obviously traceless. We are 
particularly interested in the time components which can be determined as 

Anp = 47rT;' = ^0*+l'+l"^'-n+2|$ir 
= -An{T; + T^ + T^) 

4-7^; = -^o^+i'-i^-n> (17) 

where p is the energy density. The 6 and (f components are not important for our 
analysis. Noting that the necessary non-vanishing Christoffel symbols are FJJ^ = 
= Tg^ = = a I a where the overdot represents derivative with respect to 
r], we may write the conservation equation as 9^ (y^T^) = —di (^^/gT^), where 

g = — det ga/3 = sin 6* j = a^g^ where g^ is the determinant of the 3-space; 
the terms with the F's are not written as their sum vanish due to the tracelessness. 
Integrating over the three space we find 

d, (Ea^) = -a^J d, (V^T;) = -a^J rfS.T; = J {\^+\' - |$„|2) , 

(18) 

where E is the total energy contained within a spherical volume of radius r, and 

2 

dUr = j3 sin OdOdif is an elemental area normal to the radial direction. Using the in- 
version relation, Eq. (ini), weget A2 _ = [\D_Zf\Y+f-\D+Zf\Y_f] 



dR 



\uj^\R\^ 



{\Y^\-\Y-n 



dr, 

+iuW[R*,R] + where W is the Wronskian. As the F' s are nor- 

malized to unity, when we perform the angular integration, the first term on the left 
hand side will vanish. Now R and R* satisfy the same equation and are proportional 
to the same function so that the Wronskian will also vanish. Hence Ea'^ remains con- 
stant, and so does pa^, during the evolution of the Universe. As a consequence, 
there is no net fiux of electromagnetic energy fiowing out of a closed surface. This is 
a direct proof of the well known result that the energy density of massless radiation 
scales as a~^. 
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3 Dirac Field 



In NP formalism, the four components of the Dirac spinor field of mass M are 
represented by the four functions P^,P^, — Q° and that satisfy the following 
Dirac equations [TT]: 

{D + e-p)P' + i5* + n-a)P' = 

(A + /.-7)P^ + (5 + /3-r)P° = -^Q°', 



{D + e*-p*)Q'' + {5 + 7c*-a*)Q'' = -!^pi, and 

v2 

(A + /i*-7*)Q^' + (5* + /5*-r*)Q°' = ^P°. (19) 



Substituting (^po, '-fP\ -^QO', ^Qi') with (<l>_i/2F-i/2, $1/21^1/2, $-1/2^/2, $1/21^-1/2) , 
using the directional derivatives and spin-coefficients given in Sec. 1, and following 
similar techniques, the angular parts of Eqs. (fT9|l are found to separate into 

(^L± + ^ cot 9^ y±i/2 = TAr^i/2, (20) 

while the radial and time parts satisfy 

rD±<l>±i/2 = (/A T iMra) $^1/2. (21) 

Substituting for F+1/2 = {L^ + |cot^^) F_i/2 from Eq. (f2n| into the one for y_i/2, 
and a similar process with F_i/2, gives us the two angular equations 

1 d . ^ d , icosO d 1 9^ 1 1 ^^A,, „ ,x 

■^"^^7^ ± —2^7^ + -2^7^ - -cot^^^ - - + AM = (22) 



sin 6 06 36 sin 6 dcp sin 6^ (9^?^ 4 

where the y's are the appropriate spin weighted spherical harmonics, and the eigen- 
value can be identified as A = (/ + 1/2); / is again the total angular momentum 
including the spin, i.e., / = n + 1/2, where n = 0,1, 2, 3.... The radial and time parts 
of Eq. (f2T|l can be separated with the substitution Z± = $+1/2 ± $-1/2 = R±T± to 
give 

d J\ ^ I f d 



flA*.^'^j^* = ^iU±"'V^-=*' '''' 

whence the separation constant k = pa may be identified with the co-moving mo- 
mentum. The coupled first order equations in Eq. (f23|l can be combined into pairs 
of second order equations 

± A^) T \-] R±=(-^± tMa] (-^ T iMa] T± = -k^Z^. (24) 



7 



3.1 Radial Solutions 



The radial parts take the form of the barrier penetration problem 



drl ~^ 



R± = V^Ri, (25) 



V^ = X'^±X^(^). (26) 



where the potentials are 



dr^ \ r / 

Explicitly, the potentials in the three cases are 

r (/ + l/2Tl)/r2, K = 

V± = {l + l/2)xl (/ + l/2Tcosr,)/sin2n, K = l . (27) 
[ (/ + 1/2 =F cosh r*)/sinhV^, K = -1 

In contrast to the case of the electromagnetic field, the two spin states of the Dirac 
field see different potentials. These potentials are shown in Fig. El In the flat FRW 
universe, the potential barriers basically go down as but V+ is seen to vanish 
for / = 1/2; this lowest spin state behaves as a free field as far as the radial co- 
ordinate is concerned. In the closed universe, the potentials are again infinite wells; 
for / = 1/2 however, V+ flattens out towards the origin and V_ towards vr, reaching 
the limiting value of 1/2 at the respective ends. The potential in the open universe is 
even more interesting; V+ is negative throughout for / = 1/2, reaching the minimum 
of —1/2 at the origin and approaching from below for large r,,; even for large I, V"+ 
goes down rapidly from oo at the origin, becomes negative from coshr* > (Z + 1/2), 

reaches the minimum of V+min = -|(/ + 1/2) + 1/2 - ^/{l - l/2)(/ + 3/2)j at 

coshr^, = (/ + 1/2)+ a/(/ — 1/2)(Z + 3/2) , and then rises up to zero from below 
as r^, oo; the Dirac fields of different / in the open universe see the attractive 
potential with location of the minimum at different r^, for different /. The radial 
equations turn out to be solvable for all the three cases {K = 0, ±1), and the 
behaviour of the radial part is the same for both the massive and massless fields. In 
the flat universe with K = 0, Eq. (f25|l has the solutions 

R± = r [Ay2jiTi/2{kr) + Bi/2yiTi/2{kr)] . (28) 

For regular solution at the origin, we would choose the j's, some modes of which 
are illustrated in Fig. |31 In the closed Universe, the radial equations are easily 
recognized as that for the Jacobi pj;"''^^ functions. Thus, the solutions that are 
regular between < r^, < vr can most conveniently be written as 

= Cv^sin'r. (cot ^)^'^VSi/r^/^''-^^/^^^/^)(cosr.). (29) 

Some typical behaviour are shown in Fig. |31 The solution in the open universe that 
are just the replacements k ^ ik and ir\ in Eq. (|29|l . are also shown in Fig. 

1 
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3.2 Time Dependence 



The time dependence of the Dirac field is determined by the equation 




(30) 



which corresponds to a classical oscillator in a complex time potential 



M^a^ ± iMa. 



(31) 



The mass that has coupled with the scale factor becomes dominating at large a. 
The behaviour of the solution is controlled by the evolution of the scale factor, viz., 
the Friedmann equation. Friedmann equation for a universe filled with radiation, 
massive particles and also posessing a cosmological constant can be written as 

= VL^a^ + (1 - VLm -^r- ^k) a" + VLm^ {I - v'') a'' + v'' + Qr, (32) 

where an overall factor of on the right hand side has been set equal to unity 
so that r] is in units of Hubble time; flR,M,A are the respective radiation, matter 
and vacuum densities in units of the critical density pc, and v is the rms speed 
of the matter particles, all at a reference conformal time rjo when the scale factor 
a = = 1; as before, the overdot denotes derivative with respect to rj. In this paper, 
we will concern ourselves with Qa = 0. In the form in which we have written the 
Friedmann equation, we may as well set Qr = and consider only Qm', then v = 1 
corresponds to radiation, f = is dusty matter of zero pressure, and < v < 1 
for any real matter of finite mass; this is how Vr are displayed in Fig. El for the 
three cases of f2o = 1 (flat), f2o > 1 (closed), and JIq < 1 (open) universes. For the 
massless field, Vr = 0, and the time dependence is the usual e~'^''^. In terms of the 
time marker a rather than r], Eq. ( IHUll reads as 



which is in a more convenient form to tackle the massive case. Given the time 
evolution of a by Eq. (1321), we can find series or numerical solution of Eq. (p^Hll for 
any general case. 

Here, we present the solution, in terms of Whittaker function, in the radiation 
filled flat universe where d = 1. In contrast, Sharif |H1 had to solve for the coefficients 
of a series solution. This may be due to the choice of different bases; the time 
dependence, Eq. ((201), slightly different than that of Sharif and ZeccajZj. The 
solution is 



where U and W are the confluent hypergeometric and Whittaker[T3] functions re- 
spectively. The logarithmic solution is chosen because it reverts to the massless 



d^— + -r + k^ + M^a^ T iMa T± = 0, 



(33) 




(34) 
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behaviour e~^''^ in the large k limit; viz., T+ — > A+i/vre"*'^^'^^"^**^ as — > oo. Some 
typical cases of the solution are plotted in Fig. IHl 

At early times, for any kind of FRW universe, we may make an expansion of the 
scale factor as a — > d/r/, where d/ is the slope of a at the big bang. From Eq. p2|) . 
we see that d/ ~ y/^M^ + ^r, and a is linear in rj. So the above solution, Eq. pH) . 
is valid with the replacement M MyJVtMV + VLr. 

One can factor out a time dependence T± = e=''*^'^*~*"^Ti± to write Eq. (jHUll in 
the form of a damped oscillator, 

^ ± 2^Ma^ + k'T,^ = 0, (35) 

with a velocity dependent imaginary damping that is itself time dependent. A 
formal series solution of this equation, in terms of x = rj — rjo , may be written down 

oo 

by firstly making a Taylor expansion of the scale factor as a (r/) = J2 ^^") where 

n=0 

o-n = j^n^ \rio- As di] = dx, Eq. (jT^Kll looks the same with rj ^ x. Then the series, 

oo 

Ti = ^6„x", can be substituted in the equation to find the recursion relation 

n=0 

n-2 



n{n - l)hn = {-tkfhr.-2 - {±2iM)Y^— (n - m - 1) 6„_„_i; (36) 

^-^ ml 

m=0 

Ti± will be appropriately normalized, according to the requirements that will be 
discussed below, with the choices b±o = ±l/\/2 and b±i = ^ik/\/2. For M = 0, the 
solution is Ti = e~*''^7v^. 

3.3 Dirac Particle Current 

With the generalized form of the Pauli matrix in terms of the null-tetrad 

1 / 



seen 



and the particle current written as ^ J'^ = cr^^, ^P^P^' + Q^Q^'^ , it is easily 

that the continuity equation (y/gJ^) = has to be satisfied. In particular, we 
can calculate 

J" = T^A\Z^\" + \Z-f]{\Y^f + \Y-n and 



r = -J^AZIZ^ + Z^Z1]{\Y^\' + \Y^\'). (37) 



4r2a4 

The time component of the particle current can be identified with the number density 
as n = J* = aJ^, so that na^ = a^.P is the particle number density per unit co- 
moving volume. As a consequence of the Dirac equations, Eq. (f23|) . we have on 
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hand a number of Wronskian like conditions on the radial and time parts, of which, 
two useful ones are 

d(R^R-) d(R-i-R*) r, ,9 , ,91 

^ = ^ + =ik \\R+f-\R_f] and 

dr^ dr^ 

1 \rTi |2 7 \rri |2 

= -^-^ = ik \T1T^ - T+Tl] . (38) 
dr] dr] 

From the first we see that R\R- + R+R*_ = Ci is a real constant, which must be 
equal to zero as will be shown presently; hence R\R- must be purely imaginary. 
Similarly, from the second one we find that |T+|^ + |T_|^ = C2 must be a positive real 
constant which can be normalized to unity. Now, 9,, + |^-|^] = |-R+|^ ^^^^ + 

''"'i'^' dfJ' - ^ = i'^'^- Thus. a,^j- = ^,d,U-) 

4lkf '^^^dr^ ^ '^^A?^ + l^-l^), and the particle current out of a co-moving spatial 

volume of radius r is the integration over the spatial volume: 



dr] Jq Aikf dr^, dr] 2ik dr] 

(39) 

here, we have used the property that the product of the radial functions vanishes 
at the origin, and N is the total number of particles remaining within the vol- 
ume. Thus, the net change in the number of of particles within a co-moving vol- 
ume of radius r between the times r]Q and r] is given by A^v = Nqo^q — Na^ = 
R+R- rjji^ {r]Q)f — |T_|. (77)!^]. The co-moving particle current of Eq. (f39|l has to be 



2ik 

equal to the particle flux out of a closed surface enclosing the above volume. This flux 
can be calculated by noting that Z*^Z_ + Z+Z*_ = [R\R_TIT_ + R+R*^T+T*] = 

[RIR_T*T_ + {Ci - T+Tl] = '^^If + CiT+T*. Hence, the flux is 

— J dJ^rJ^ = ^2ik~ '^^'^ri^ ~^ (term with Ci). For this flux to be equal to the particle 
current of Eq. as required by continuity equation, Ci = 0, as was asserted 

earlier. Thus the particle current out of a co-moving volume is finite, unless \T^f 
is a constant, as would be for the massless case where the time dependence e^*^'' 
keeps the co-moving particle number conserved for all times. 

For the massive case on the other hand, as |T+| cannot be a constant, the particle 
current does not vanish, but its magnitude goes down nonetheless with the expan- 
sion. Explicitly, in the radiation filled flat case that was solved in Eq. the ap- 



propriately normalized time functions are T+ = J ,2\'^^ / ,2 \ e ( tit + 

y cosh(ff)+sinh(^) VM 

and T_ = / / ,2n^ ^^e~^^'''i'' I'^U iMr]A whose sum of the modulus- 

\/ cosh(#Usmh(#) \iM^2^ I ) 



squares is unity, while R\R- = i\A+f r'^ji_i/2{kr) ji+i/2{kr); then the particle 
currents can be easily evaluated. A point to be noted is that in the flat universe, 
taking the limit r — > cxd of Eq. (jSHl), the product of the radial functions vanishes, 
so that the co-moving particle number of the whole universe still remains constant. 
The behaviour of dr, (Na^) are shown in Fig. Elas a function of r at a time r] = 1. 



11 



As is evident from the relation itself and from the figure, for a given k, the current 
becomes zero at the zeroes of the j's at particular values of r; these nodal points 
are not dependent on the mass of the Dirac field, but only on k and /. Hence, the 
current oscillates between positive and negative values. This is also seen in Fig. 

As a result, the particles flow out of regions of negative current into those of 
positive, and the radial dimensions develop areas of depleted and excess particles. 
The magnitude of the current at a particular r just goes down with rj, Fig. ((Hj). It 
is seen that the magnitude of the current also goes down with the mass. 

This flow of electrons and protons in or out of different regions of the universe 
would constitute a current. Such a current would be very prominent during the time 
between the electron-positron annihilation which left an excess of electrons and the 
recombination when the electrons were absorbed in the neutral atoms. It is reason- 
able to expect that such charged current must have generated a cosmic magnetic 
field. The magnetic field would greatly modify the cosmic hydrodynamics, and have 
immense consequence on distribution of matter and structure formation. Besides 
the charge current, massive neutrinos would also similarly flow between different re- 
gions of the Universe. Consequently, homogeneity would be destroyed, and regions 
of excess and reduced number of neutrino would appear. Such inhomogeneity in 
neutrino distribution could also greatly affect the process of structure formation. 

3.4 Energy-Momentum of Dirac Field 

We now look into the behaviour of the energy momentum tensor. The energy- 
momentum spinor given in Ref. can be used to evaluate it in tensorial form; 
the components that we are particularly interested in are: 

• f2 

= -Tr=^ [^l/2D-n/2 - $-l/2/^+$-l/2 - CC] (iFl/^r + |>^-l/2r) , 

Te' = -T^. (40) 

The angular components are not explicitly required for our purposes; it suffices us 
to know that Tl^ is traceless. The conservation equation then reads drj [^/gTJ|) = 
—di [y/gT^), and is just the energy density p. The angular integration may be 
easily done. Evaluating the square bracket in the expression for we find that it is 
—ik (|T+|^ + \T^f) [R^R- + R+R^]; it was shown in the previous subsection that 
the time term in the round bracket is constant while the radial term in the square 
bracket vanishes. Hence the momentum fiux out of a co-moving volume is zero and 
a^p remains constant during the evolution of the Universe. It is indeed intriguing 
that the Dirac particle current is finite and the number density is not conserved 
within a co-moving volume, but the momentum current vanishes and the co-moving 
energy density remains constant. 
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4 Conclusions 



In this paper, the relevant fields of electrodynamics, viz., the free Maxwell and Dirac 
fields, in Friedmann-Robertson- Walker spacetime have been investigated using the 
Newman-Penrose method. All the variables are found to be separable, and the 
angular solutions are the spin-weighted spherical harmonics. The massless fields 
have the usual exponential time dependencies. All the radial parts reduce to the 
quantum mechanical barrier penetration problem, with well behaved potentials that 
are basically the centrifugal energies. The potentials seen by one component of the 
Dirac field, are interesting; its lowest angular momentum state sees no potential 
in the flat universe, while it sees an attractive one throughout the open universe; 
from afar, all angular momentum states of this component see attractive potential 
in the open universe. Consequences of this effect may provide a means to determine 
whether the Universe is flat, open or closed. All the radial equations are solved. 

The time evolution of the massive Dirac field is found to be coupled to the 
evolution of the cosmic background. Although the temporal part of the massless 
Dirac field is that of a free classical oscillator, the massive field experiences a complex 
potential; the real part of this time-potential is a downturned parabola as a function 
of the scale factor, going to deeper negative values with the expansion; the imaginary 
part is proportional to the rate of change of the scale factor. However, in the special 
case of the radiation filled flat universe where a = 1, an analytic solution has been 
determined. 

The behaviour of the respective conserved currents and fluxes have also been 
discussed. The conservation of the energy-momentum of the Maxwell held leads 
to expected results: the momentum flux out of a closed co-moving surface times 

vanishes; consequently, the energy density scales as a^^. The behaviour of the 
massless Dirac held is also similar: its co-moving number density is conserved and 
there is no particle current out of a co-moving volume; its energy density also scales 
as a^^. The massive Dirac field on the other hand shows a very different behaviour: 
the co-moving particle current is not zero, and obviously, the particle number within 
a co-moving volume is not conserved; this current depends on the mass of the field. 
Hence, regions that are depleted of particles of certain angular and linear momenta, 
and others with excess, will appear in the universe. This could lead to effects that 
have not yet been considered in studying structure formation. Similar currents of 
electrons and protons that existed before the recombination era should have gener- 
ated a cosmic magnetic field. The existence and shape of such a magnetic field would 
have strong influence on the hydrodynamics of cosmic matter, resulting in hitherto 
unaccounted effects on the structure formation processes. Further consequences and 
interaction between the flelds will be investigated in subsequent works. 
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Figure 1: The tortoise coordinate function of r. In flat space with K = 0, 

(sohd), the two co-ordinates are identical. For K = 1 (dotted), r^, plateaus off to a 
value of 71 at large r. For K = —1 (dashed), the maximum value that r can have is 
2, where n oo. For small r, all three cases are alike. 




Figure 2: The potential seen by electromagnetic waves. / = 1, 2, and 3, for the solid, 
dotted and dashed lines, i^' = 0, 1, and —1 from left to right. The potential is a 
wall in the flrst, become inflnitely high wall at each boundary in the second, and an 
exponentially decaying wall in the third case. The walls are thicker for higher I. 
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Figure 3: The potentials seen by the radial part of the Dirac field. / = 1/2,3/2 
and 5/2 for the solid, dotted and dashed curves. K = 0,1 and —1 from top to 
bottom. The figures on the left are V+ and the right are V-. For / = 1/2,1 + 
vanishes in the flat universe and becomes a negative well in the open universe. At 
larger r* in the open case, V+ falls down rapidly from oo at r* = to a negative 
minimum value and rises asymptotically to zero from below. In the closed universe, 
for / = 1/2, 1/+ ^ 1/2 as 0, and V_ 1/2 as n. 
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Figure 4: The radial parts of the Dirac field. K = 0,1 and —1 from top to bottom. 
In the top one: k = 1/2, 3/2 and 5/2 from left to right: / = 1/2, 3/2 and 9/2 for the 
solid dotted and dashed curves. In the middle and bottom rows, / = 1/2,3/2 and 
9/2 from left to right; k = I + 1,1 + 2 and Z + 4 for the solid, dotted and dashed 
lines. In the bottom row, the minimum of are —0.5, —0.27 and —0.25 from left to 
right; in the left, k = 0.1,0.4 and 0.6 for the solid, dotted and dashed lines; in the 
other two, k — 0.1, 0.24 and 0.3 for the solid, dotted and dashed lines. 



17 




Figure 5: The potentials seen by the temporal part of the Dirac field as a function 
of the scale factor a. The rest mass has been set to unity. K — 0,1 and —1 from 
top to bottom. The real part of the potentials (solid lines) is the same. The broken 
lines are ImVr+ = — ImV^-. The dotted lines are for radiation filled universe with 
particles of rms speed vq = c, the dashed lines with vq = 0.3c are for the universe 
near transition to the non-relativistic phase, and the dot-dashed lines are for dust 
filled universe with vq = 0. 
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Figure 6: The time part of the massive Dirac field in a radiation filled flat universe. 
M = 1 in the solid and broken curves, for which k = 1 and 2 respectively. M = 2 
in the dot-dashed and the dashed curves, for which k — 1 and 2 respectively. The 
functions are normalized according to iT+j^ + \T_f — 1. 




Figure 7: The Dirac current in a radiation filled flat universe as a function of r. 
They are all shown at time rj — 1. All the curves are for I — 1/2, except the dotted 
one of 5/2. The mass is 1 for all except the dashed-dotted curve which is for M = 2, 
and k = 1 for all except the broken line for which it is 2. The current is seen to go 
down with M. The particles flow out of the regions of negative current, and into 
those of positive. So the regions develop into those depleted of certain k and /, and 
those with excess. For example, the region up to r = 3 is almost devoid of particles 
of Z > 1/2 



19 



a^[Na3] 
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Figure 8: The time dependence of the Dirac current in a radiation filled flat universe 
at particular r. All are at r = 4, except the dashed-dotted one at r = 1. Values of 
/, M and k are as in the Fig. ^ . The magnitude of the current goes down with 
time. 
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